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We present Gemini optical spectroscopy of 23 young star clusters in NGC 3256. We find that 
the cluster ages range are from few Myr to ^ 150 Myr. All these clusters are relatively massive 
(2-40) xlO^ A^0 and appear to be of roughly 1.5 ZQmetallicity. The majority of the clusters 
in our sample follow the same rotation curve as the gas and hence were presumably formed in 
the molecular-gas disk. However, a western subsample of five clusters has velocities that deviate 
significantly from the gas rotation curve. These clusters may either belong to the second spiral 
galaxy of the merger or may have formed in tidal-tail gas falling back into the system. We discuss 
our observations in light of other known cluster populations in merging galaxies, and suggest that 
NGC 3256 is similar to Arp 220, and hence may become an Ultra-luminous Infrared Galaxy as 
the merger progresses and the star-formation rate increases. 

Some of the clusters which appeared as isolated in our ground-based images are clearly resolved 
into multiple sub-components in the HST-ACS images. The same effect has been observed in the 
Antennae galaxies, showing that clusters are often not formed in isolation, but instead tend to 
form in larger groups or cluster complexes 
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ABSTRACT 



In the prevailing picture of hierarchical early- 
type-galaxy formation, small fragments form first 
and then merge into larger and larger pieces until 
the system resembles a large, smooth, anisotrop- 



globular clusters has an important role to play in 
testing the predictions of this theory and in an- 
swering the question about if and when the bulk 
of this merging took place. Studies of the Galactic 
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globular cluster system have been fundamental to 
the development of ideas on how the Galactic halo, 
and perhaps the entire galaxy, has been assem- 
bled from merging fragments (Searle & Zinn 1978; 
Da Costa & Armandroff 1995). Studies of extra- 
galactic cluster systems have revealed evidence for 
multiple populations of clusters and allow us to 
use them to connect the various phases of galaxy 
evolution. 

GCs in galaxies more distant than about a 
megaparsec cannot be resolved into stars, even 
with HST. However, the fact that they contain 
simple stellar populations means that their inte- 
grated properties can be used to determine mctal- 
licities and ages. Surveys of large numbers of glob- 
ular cluster systems reveals that >90% of all Es 
have bimodal color distributions, and only 10-20% 
in SOs (Kundu & Whitmorc (2001a, 2001b)). Al- 
most all galaxies host a GC population with a peak 
in their color distribution at (V-I)=0.95. Avail- 
able metallicity measurements suggest that these 
are old, metal-poor clusters like the Galactic halo 
GCs (Geisler, Lee, & Kim 1996). The second peak 
is redder than the first and implies that the red 
GCs are more metal-rich. In general, these GCs 
are also old, but the age-metallicity degeneracy in 
broad-band colors makes age differences hard to 
determine unless the GCs are rather young (see 
Whitmore et al. 1997). 

There are several competing scenarios for how 
the globular cluster sy stems of galaxies formed. 



Ashman fc Zepf I (;1992) predicted that the merger 
of two disk galaxies would produce an elliptical 
with a bimodal GC color distribution. The blue 
clusters are Galactic-halo like clusters from the 
progenitor galaxies. These clusters would have 
formed like those in the Galaxy, perhaps from the 
accretion of dwarf galaxies with large numbers of 
GCs. Dwarf elliptical galaxies were efficient at 
forming clusters and these clusters resemble those 
in the Galactic halo (Miller et al. 1998). In the 
Ashman & Zepf scenario, the red population is 
formed during the merger process from the metal- 
enriched gas in the disks. The formation of new 
clusters also alleviates the problem that ellipticals 
have specific globular cluster frequencies (number 
per unit luminosity) about a factor of two higher 
than spirals do (Schweizer 1987). An alternative 
scenario is that all the clusters were formed "in 
situ," in a multi-phase collapse of a single poten- 



tial well (Forbes et al. 1997) similar to the early 
monolithic collapse picture of the formation of the 
Galaxy (Eggen, Lynden-Bell, & Sandage 1962). In 
this case, the metal-poor clusters form first in the 
halo and and metal-rich clusters form later dur- 
ing the final collapse, or recoUapse, stages from 
metal-enriched gas. Both these scenarios fit into 
the overall picture of hierarchical galaxy forma- 
tion, but they differ in when the bulk of the merg- 
ing takes place. 

An important result of the merger scenario is 
that the merger of two spiral galaxies can cause 
the formation of many young star clusters in the 
starburst resulting from the collision of the two gas 
disks (Schweizer 1987; Ashman & Zepf 1992). The 
creation of star clusters may alleviate the prob- 
lem that elliptical galaxies have specific globular 
cluster frequencies {Sn, the luminosity- weighted 
number of GCs) about a factor of two higher than 
in spirals. The evidence is growing that signifi- 
cant numbers of star clusters are formed during 
galaxy mergers. Hubble Space Telescope observa- 
tions of NGC1275 (Holtzman et al. 1992; Carl- 
son et al. 1998) were some of the first to find 
very luminous blue objects with the properties of 
young GCs in a galaxy that may have had a re- 
cent merger. Large young cluster populations with 
ages between 10 and 500 Myr were soon found 
in obvious merger remnants like NGC7252 (Whit- 
more et al. 1993; Miller et al. 1997), NGC3921 
(Schw eizer et al. 1996), NGC3256 (|Zepf et al.l 
19991) . NGC4038/4039 (Whitmore & Schweizer 
1995; Whitmore et al. 1999). Spectra of the 
brightest young clusters in NGC1275 (Zepf et 
al. 1995) and NGC7252 (Schweizer & Seitzer 
1993,1998) confirmed that the clusters were be- 
tween 0.5 and 1 Gyr old with roughly solar metal- 
licities (Figure 2). Since these clusters are many 
internal crossing times old, they would seem likely 
to survive to become the red GC populations of 
the faded merger remnants (e.g., Whitmore et al. 
1997; Goudfrooij et al. 2001a). 

These observations allow us to proceed to the 
next level of understanding the evolution of glob- 
ular cluster systems; we must understand in more 
detail how star clusters form and how systems of 
star clusters evolve. It is thought that star clusters 
are formed in high-pressure environments (Harris 
& Pudritz 1994; Elmegreen & Efrimov 1997), be 
they caused by the collisions of giant molecular 
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clouds or by a general pressure increase in the 
gas surrounding molecular clouds. Observations 
of the Antennae can now provide some key pa- 
rameters about the state of the ISM during clus- 
ters formation and the feedback produced by the 
young clusters. Zhang, Fall, & Whitmore (2001) 
compared the locations of different cluster popu- 
lations with emission from the ISM at wavelength 
from X-rays to radio and found that the youngest 
clusters are associated with molecular cloud com- 
plexes and may lie in regions of high HI velocity 
dispersion. However, the small velocity dispersion 
of the clusters among themselves strongly suggest 
that it is not high-velocity cloud-cloud collisions 
that drive cluster formation, but the general pres- 
sure increase experienced by gas during the merger 
(Whitmore et al. 2005). Feedback, seen in the 
form of Ha bubbles around young-cluster com- 
plexes, may enhance this process. 

An interesting question is how globular clus- 
ter systems evolve. Most very young star-cluster 
systems studied in merging galaxies have power- 
law luminosity functions of the form M~^ (e.g. 
Schweizer et al. 1996; Miller et al. 1997; Whit- 
more et al. 1999; Zhang & Fall 1999). However, 
the luminosity (or mass) function of old globular 
clusters has a log-normal shape with a peak at 

7.3 (^ 105 M©; Harris 1991). Therefore, 

either the initial mass function of star clusters was 
different in the past, perhaps due to low metallic- 
ity, or a substantial fraction of the young clus- 
ters must be destroyed for the young mass func- 
tions to evolve into what we see in older systems. 
A great deal of theoretical work has gone into 
globular cluster destruction processes (e.g. Fall 
& Rees 1977; Gnedin & Ostriker 1997; Vesperini 
1997/98; Fall & Zhang 2001). The main processes 
that can destroy globular clusters are stellar evo- 
lution, 2-body relaxation, dynamical friction, and 
disk shocking, or if the clusters formed in gas rich 
environments, interactions with GMCs can play a 
significant role (Gieles et al. 2006). The models 
of Fall & Zhang (2001) predict that there may be 
radial variations in the peak of the mass function 
within a galaxy and that the peak will shift to 
higher masses with time. Many of these processes 
depend on the relative velocities of the clusters 
and field stars and on the cluster orbits. 

In this paper we present some first results of 
a large spectroscopic survey of star clusters in 



merging and interacting galaxies. We focus on 23 
bright star-cluster candidates in the main body of 
NGC 3256 observed with GMOS on Gemini South. 
Other targets in our survey, to be presented in fu- 
ture papers, include NGC 4038/39 (The Anten- 
nae), NGC 6872, Stephan's Quintet, and M82. 

NGC 3256 was classified as an intermediate- 
stage mer ger in Toomr e's list of nearby merging 
systems (jToomre II1977I ). The merger is more ad- 
vanced than systems like NGC 4038/39, in which 
the two original disks are still distinct, but the two 
nuclei (Moorwood & Oliva 1994; Norris & Forbes 
1995; Lira et al. 2002; Neff et al. 2004) have not 
merged yet either. They have a separation in pro- 
jection by 5" , or - 900 pc. Hence, NGC 3256 is 
not as relaxed a system as NGC 7252 is. The outer 
parts of the system are characterized by shell-like 
features and two extended tidal tails that are typ- 
ical of merging galaxies. The body of the system 
is criss-crossed by dust lanes that enshroud an on- 
going starburst: the far-infared luminosity. X-ray 
luminosity, and star formation rate are the highest 
of all the systems in the Toomre sequence. This 
starburst has c reated over 1000 star clusters in the 
central region (IZepf et al.lll999l) as well as in the 
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tidal tails (jKnierman et al 
i2007h . 

The current paper is organized as follows: Sec- 
tion [2] describes the observations. In § [3] and § |4] 
we derive the ages, masses, extinctions, metallici- 
ties, and line-of-sight velocities of the 23 clusters. 
Finally, § [5] discusses and summarizes the results. 

NGC 3256 is located at ajaooo = 10'^27'"51!3, 
^.nooo — — 43°54'14" and has a recession velocity 
relative to the Local Group of cz^^^.^ = -1-2804 ± 6 
km , which places it at a distance of 36.1 Mpc 
for Hq — 70 km s~^ Mpc~^. At that distance, 
adopted throughout the present paper, 1" = 175 
pc. The corresponding true distance modulus is 
(to— M)o — 32.79. Because of the low galactic lati- 
tude of NGC 3256, 5 = +11?7, the Milky Way fore- 
ground extinction is relatively high. Ay — 0.403 
(Schlegel et al. 1998), whence the apparent visual 
distance modulus is — My = 33.19. 

2. Observations and Reductions 

Imaging and spectroscopic observations of star 
clusters in NGC 3256 were made with GMOS-S 
in semesters 2003A and 2004A. The data were 
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obtained as part of two Director's Discretionary 
Time programs, GS-2003A-DD-1 and GS-2004A- 
DD-3. Our images cover the typical GMOS-S field, 
which measures approximately 5.'5 x 5f5. They 
were obtained through the g' and r' filters. Four 
CMOS masks with slitlets were used for the spec- 
troscopic observations. We used the B600 grat- 
ing and a slitlet width of 0'.'75, resulting in an 
instrumental resolution of 110 km/s at 5100 A. 
The spectroscopic observations were obtained as 8 
individual exposures with exposure times of 3600 
sec each. Our spectroscopy of 70 cluster candi- 
dates yielded only 26 objects that were bonafied 
star clusters in NGC 3256. Of these, three are 
located in the western tidal tail and have been de- 
scribed in Trancho et al. (2007). In the present 
paper we focus on the 23 star clusters located in 
or near the main body of NGC 3256. Tabic [T] lists 
these star clusters. Column (1) gives the adopted 
cluster ID, columns (3)-(4) the coordinates, and 
columns (5) and (6) the absolute magnitudes Mg' 
and Mr> and their errors. The magnitudes have 
been corrected for Galactic extinction, but not for 
any internal extinction. 

Figure [T] shows an HST/ AGS image of the main 
body of NGC 3256, with the observed candidate 
clusters marked by their ID numbers. 

The basic reductions of the data were done us- 
ing a combination of the Gemini IRAF package 
and custom reduction techniques, as described in 
Appendix A. 

3. Derivation of Cluster Properties 



The derivation of cluster properties (such as age 
and metallicity) based on the strengths of stellar 
absorption lines through optical spectroscopy is 
not straightforward, due to degeneracies between 
age, metallicity, and extinction. Multiple stud- 
ies have addressed this problem (e.g., Schweizer & 
Seitzer 1998; Schweizer, Seitzer, & Brodie 2004; 
Puzia et al. 2005), and here we extend previous 
studies. 

Although in some of our cluster spectra the 
strengths of stellar absorption lines cannot be 
measured due to strong emission lines, the fluxes/equivalents 
width(EW) of the emission lines of the surround- 
ing H II region can be measured. In these cases 
a chemical abundance can be estimated for the 
H II region, in which the cluster has recently 



formed, from line-emission measurements (e.g., 
Kobulnicky & Kewley 2004; Kobulnicky & Phillips 
2003; Vacca & Conti 1992). This abundance is ex- 
pected to be the same as that of the young stellar 
cluster itself, and hence complements abundance 
measurements of the absorption-line dominated 
clusters. 

Below, we outline the method adopted in the 
present study. 

3.1. Extinction, Age, and Metallicity 

3.1.1. Absorption-Line Clusters 

Wc first select the model spectra by Bruzual & 
Chariot (2003, hereafter BC03) and by Gonzalez- 
Delgado et al. (2005, hereafter GD05) . Then we 
smooth the model spectra to the same resolution 
as the observed spectra. Then we compare the 
cluster spectra with the models for clusters of so- 
lar metallicity, to which we have applied various 
amounts of extinction (using de Galactic extinc- 
tion law for Savage & Mathis 1979) , Ay = 0-10 in 
steps of 0.1 mag. We select the best fitting model 
via minimized x^, Modelbcst (age,>lv'), and correct 
the observed spectrum for the derived extinction 
Ay to yield Clusterobs.oxt- 

This procedure was carried out for the BC03 
and GD05 models independently, and we note that 
for individual clusters the results are very similar. 
Due to the finer grid of young ages in the GD05 
cluster models we adopted these for our further 
analysis. 

We then inserted the extinction-corrected spec- 
tra into the IDL implementation of the Penalized 
Pixel Fittingroutine (pPxF) of Cappellari & Em- 
sellem(2OO40 This routine determines the best 
linear combination of template spectra plus an an- 
alytic polynomial to reproduce the observed spec- 
tra and returns, in addition, the radial velocity 
of each cluster. For template spectra we used 
all available models (both in age and metallic- 
ity) of GD05. We suppressed the use of any ad- 
ditional polynomial in order to preserve the con- 
tinuum shape of each observed cluster. Emission 
lines and regions of the spectra affected by ar- 



^We used the GD05 spectra which have a resolution of 0.3 
A at 5100A, corresponds to cr 0.3/5100 * 3e5/2.35 = 7.5 
km/s and our spectra have a cr ~ 3.96/5100 * 3e5/2.35 = 
99.12 km/s. The quadratic difference is sigma necessary to 
use in ppxf is cr = ^{99.12^ - 7.5^) = 98.83 km/s 
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tifacts were masked during the fits, though care 
was taken to minimize the number of such areas. 
Via this procedure we obtained a template spec- 
trum for each of our clusters, Clustertcmp, as well 
as an emission spectrum (the difference between 
Clusterobs,oxt and Clustertcmp)- 

In order to find the age and metallicity of each 
cluster, we measured the line strengths of the 
hydrogen Balmer lines as well as of prominent 
metal lines (See Table 2). For this we used the 
Lick indices (Faber et al. 1985; Gonzalez 1993; 
Trager et al. 1998) as well as the indices defined 
by Schweizer & Seitzer (1998) for young stellar 
populations. To make the measurements, we used 
the routine Indexf (Cardiel et al. 1998), which 
finds the line strengths and errors by performing 
Monte-Carlo simulations on the spectra, using in- 
formation derived from the error spectra (i.e., the 
placement of the continuum bands and noise in 
the data). Indexf was run on Clustertcmp instead 
of Clusterobs,cxt- The reason for this is that we 
found that the measured index strengths of the 
Clusterobs,cxt depended on the S/N ratio of the 
spectrum. This is shown in Fig. [3l where we plot 
the nine measured indices of an observed cluster 
(T130 in The Antennae, which is our best S/N 
cluster and has identical setup and was observed 
in the same way) degraded to various S/N ratios. 
The open symbols represent the measurement of 
the line strength for each index when Indexf is run 
on Clusterobs.oxt, while the filled symbols represent 
the measurements carried out on Clustertcmp- The 
lines show the average index strength for the five 
highest S/N measurements. From this numerical 
experiment we conclude that the measurements on 
Clustertcmp reproduce those of Clusterobs,cxt for 
high S/N, but remain accurate to S/N <10, while 
measurements on Clusterobs,cxt begin to show sig- 
nificant scatter below S/N « 15. 

However, we note that with this adopted pro- 
cedure the measured H/3 index is always system- 
atically off. The models never reproduce an ab- 
sorption feature on the red side of the line. Addi- 
tionally, similar tests were performed where we in- 
serted the model spectrum (i.e. treating the mod- 
els as observations, degrading them in S/N and 
finding the indices). These tests showed a system- 
atic offset in the USa measurements although the 
other lines were well reproduced. Therefore, we 
removed both the H/? and USa indices from our 



list for further analysis. 

Using Indexf we also measured the line strengths 
of the GD05 models, for all model ages and metal- 
licities. The age and metallicity of each observed 
cluster was then determined by comparing its line 
indices to that of models, weighted by the respec- 
tive errors, in a least sense. In order to test 
the robustness of this technique we added random 
errors to the measured indices (using a normal dis- 
tribution with a dispersion corresponding to the 
l-cr measurement error) and re-did the analysis. 
This was done 5000 times for each cluster. The 
final age was then determined by creating a his- 
togram of the derived ages and fitting a gaussian 
to it (in logarithmic space), where the adopted age 
is the peak age of the gaussian and its error is the 
standard deviation. The cluster metallicity was 
found by simply averaging the derived metallicity 
of the 5000 runs. Examples of this process are 
shown in Fig. [51 

Finally, to check the consistency of our results, 
we plotted the spectra of Clusterobs.cxt and the 
best fitting model (i.e., the model closest in age 
and metallicity). If the fit was not satisfactory, 
then we began the entire process again, eliminat- 
ing the initially derived extinction from the op- 
tions. This was the case for only a handful of 
young clusters that were initially fit with low ex- 
tinctions and higher ages. Two examples of cluster 
spectra are shown in Fig. [S] 

The black lines are the Clusterobs.cxt spectra, 
while the red and green lines represent the best 
fitting model (age, metallicity) and the residual 
(Clusterobs.cxt — template — constant). T1002 in 
the right panel is clearly very young and has, as 
such, still ionized gas associated with it. In the 
observed and residual spectra, we clearly see emis- 
sion fines from H7, H/3 and [OIII]AA4959,5007. 

As a further test of the derived cluster ages, 
we chose a subset of indices that are good at dis- 
tinguishing between the age and metallicity of a 
cluster. For illustrative purposes, Fig.[4]shows the 
[MgFe] index (Thomas et al. 2003) plotted versus 
the H7 index for both the GD05 and BC03 mod- 
els. The indices from the models are also shown. 
From these diagrams we check for consistency in 
the derived ages and metallicities of the clusters. 

The derived ages, extinctions, and metallicities 
are given in Table [TJ 
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3.1.2. Emission-Line Clusters 

For the youngest clusters with Httle or no ab- 
sorption features in their spectra the task is much 
easier. First, we assign ages to these clusters of less 
than 10 Myr, due to the presence of large amounts 
of ionized gas around the cluster. Age dating can 
be refined to some degree by the presence or ab- 
sence of Wolf-Rayet features (e.g., Cluster T2005, 
see Fig[S]); however, that is beyond the scope of 
the present paper. The extinction of these clus- 
ter is calculated from the H7 to H/3 emission-line 
ratio. 

We adopt the chemical analysis method from 
Kobulnicky & Kewley 2004 (hereafter KK04) 
to determine the metallicity. We measure the 
EW ratio of the coUisionally excited [OII]A3727 
and [OIII]AA4959,5007 emission lines relative to 
the H/3 recombination line (known as R23) and 
[OIII]AA4959,5007 relative to [OII]A3727 (knows 
as O32) , along with the calibrations on KK04 
(their Fig. 7 - upper brach) and the solar abun- 
dances by Edmunds & Pagel (1984). Instead of the 
traditional flux ratio, the KK04 method uses EW 
ratios that have the advantage of being reddening 
independent. 

As can be seen in Table [1] the metallicites found 
for absorption-line and emission-line clusters agree 
well, giving us confidence in the robustness of the 
diagnostic methods and results. 

3.2. Masses 

In order to calculate the mass of each cluster 
we compared the photometry {g' and r') with the 
BC03 SSP models, assuming a Chabrier (2003) 
stellar initial mass function and solar metallicity. 
We then used the age dependent mass-to-light ra- 
tio from the models to convert our derived abso- 
lute magnitudes (observed magnitudes corrected 
for Galactic extinction, internal extinction, and 
the assumed distance modulus) to mass. Errors on 
the mass were estimated from the derived errors 
on the age and photometric errors. Systematic er- 
rors (e.g., errors associated with the distance to 
NGC 3256) are not included. Table [3] gives the 
derived masses of the clusters. 

3.3. Velocities 

In the case of absorption- line clusters, we used 
the IRAF task rvsao.xcsao for the determination 



of the redshift from the individual spectra, us- 
ing three different type (A, O, B) radial-velocity 
standard stars (HD 100953, HD 126248, and 
HD 133955) observed at the same resolution as the 
clusters. The three template stars were employed 
to reduce the systematic errors introduced by the 
effect of template mismatch when computing the 
redshift using the cross-correlation technique. 

For the emission-line clusters, velocities were 
measured from the observed emission lines using 
the IRAF task rvsao.emsao. 

In both cases the velocities were corrected to 
heliocentric (see Table 3). 

Figure [7] shows the positions of the observed 
clusters within NGC 3256 (shown in contours to 
highlight its main features), marked with the clus- 
ter metallicities, extinctions, ages, and velocities, 
respectively. 

4. Cluster Kinematics: Two Populations? 

There is strong evidence that the molecular gas 
in the central region of NGC 3256 lies in a disk 
that rotates (Sakamoto et al. 2006). The rota- 
tion axis of this gas disk lies approximately along 
the north-south direction, which is also the ap- 
parent minor axis of the main optical disk. It is 
interesting to compare the observed cluster radial 
velocities to the molecular-gas velocities at each 
cluster's position. 

Figure [5] shows the measured radial velocities 
of the clusters plotted versus the cluster right as- 
cension (RA), corresponding roughly to their pro- 
jected position along the major axis. Superim- 
posed is the rotation curve for the molecular gas 
measured by Sakamoto et al. (2006). The figure 
suggests that the majority of the clusters are still 
associated with the gas from which they formed 
(see Table [3|). At least 15 of the 23 observed clus- 
ters show clear evidence of corotating with the 
molecular-gas disk. Hence, we will refer to these 
clusters as "disk clusters." 

4.1. Origin of the Disk 

Sakamoto et al. (2006) suggest that the molecular- 
gas disk may have formed during the merger of 
the parental spiral galaxies. By using the ages of 
the disk clusters we can put a lower limit on the 
longevity of the disk. These ages range from re- 
cently formed (e.g., T761: <10 Myr) to 100 Myr 
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old (T112). Thus the molecular-gas disk must 
have existed for at least 100 Myr. 

The NGC 3256 merger probably began approx- 
imately - 500 Myr ago (EngHsh et al 2003). It 
is not yet completed, as the two nuclei have yet 
to merge (Moorwood & Oliva 1994; Norris & 
Forbes 1995; English & Freeman 2003). There- 
fore, if the observed molecular-gas disk and clus- 
ters formed during the merger, the disk must have 
begun forming early in the merger. Sakamoto 
et al. (2006) compare the NGC 3256 system to 
that of NGC 7252, a recently formed merger rem- 
nant which also hosts a molecular-gas disk. How- 
ever, in contrast to the situation in NGC 7252 the 
two nuclei of NGC 3256 have yet to merge, which 
may disrupt any current large-scale gaseous disk 
(Barnes 2002). Hence, the two merger systems 
may not presently be comparable. An alternative 
hypothesis, however, is that the observed gas disk 
belongs to one of the two original spiral galaxies, 
so that the observed disk clusters simply formed 
in that disk as part of the enhanced star-formation 
activity caused by the gravitational interaction. 

In either scenario, we would expect older star 
clusters to be present as well. The fact that they 
are not detected in the present study is readily ex- 
plained by selection bias: we selected the brightest 
clusters, which tend to be young, for spectroscopy. 

4.2. Non-Disk Clusters 

In the western section of the galaxy we find five 
clusters which have velocities apparently inconsis- 
tent with an extrapolation of the molecular-gas 
disk velocities. The clusters have ages between 
< 7 Myr (e.g. T96) and - 150 Myr (e.g. T1078). 
These clusters may belong either to the other spi- 
ral disk (which may lie behind the observed disk 
in projection) or to material which has become 
dissociated from the original disks due to the in- 
teraction/merger. These clusters are also located 
spatially near the beginning of the western tail 
(see Fig. 1 in Paper I). A detailed comparison 
with the HI position-velocity diagram of English 
et al. (2003) shows that the HI tail begins approx- 
imately 45" to the west of the observed clusters. 
However, as Fig. [TU] shows, these clusters are co- 
incident spatially and kinematically with HI gas 
that has a very different radial-velocity distribu- 
tion from that of the molecular gas. The H I radial 
velocities reach a minimum near the CO rotating 



disk and a maximum at the kinematic center of 
NGC 3256, whereupon they begin dropping again 
inside the western tidal tail. One possible inter- 
pretation is that the gas-velocity anomaly, noted 
already by Enghsh et al. (2003), is caused by gas 
falling back into the central parts of NGC 3256 
from one of the tidal tails (perhaps the eastern 
one). 

Two other clusters located closer to the ob- 
served galactic center stand out in terms of their 
kinematics. These clusters (T779 and T343) have 
velocities larger than that expected if they were 
part of the rotating molecular-gas disk (although 
T343 is only incompatible with the disk velocity 
at the 1.5cr level). Both clusters are very young 
and have only modest extinction {Ay— 0.4-0.8). 
We note that these two clusters are located in a 
part of the galaxy where there is a rather large 
scatter in the measured velocities of the clusters 
(e.g. T201, T343, T356, T779) and thus their de- 
viation may be part of a larger trend. It is possible 
that we are seeing a heating or beginning destruc- 
tion of the disk due to the interaction/merger, and 
that star-formation is proceeding from an ordered 
phase, i.e. in a disk, to that of a more chaotic 
phase where dispersion dominates over rotation. 

5. Discussion 

5.1. Environment of the Clusters 

Some of the emission line clusters that ap- 
peared isolated in our ground based images are 
clearly resolved into multiple subcomponents in 
the HST/ ACS images. The same phenomenon 
has been observed in the Antennae galaxies (e.g., 
Whitmore & Schweizer 1995; Bastian et al. 2006), 
showing that clusters are often not formed in iso- 
lation but instead tend to form in larger group- 
ings, or cluster complexes. These complexes are 
thought to be a short-lived phenomenon as they 
disperse on short timescales, although merging 
within the central parts of the groupings is pos- 
sible (e.g., Fellhauer & Kroupa 2002). The rem- 
nants of such cluster-cluster mergers are an attrac- 
tive means to form extremely large clusters, such 
as W3 and W30 in NGC 7252 (Kissler-Patig et al. 
2006). 

Hence, it is possible that in a few cases (which 
are in very crowded regions, e.g. T2005 and T116) 
we may be over-estimating the mass of an appar- 
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ent "cluster" if, in fact, it is made up of several 
clusters. 

5.2. Star/Cluster Formation Rates 

NGC 3256 contains the most molecular gas 
among the merging galaxies and merger remnants 
of the Toomre sequence (1.5 x 10^° -Mq: Casoli 
et al. 1991; Aalto et al. 1991; Mirabel et al. 1990, 
from Zepf et al. 1999). Thus, there will be plenty 
of gas left to fuel a massive starburst when the 
nuclei merge (e.g., Mihos & Hernquist 1996). At 
that time, one may expect the star/cluster for- 
mation rate to increase substantially (between 3 
and 10 times depending on the encounter param- 
eters and the time to nuclear coalescence) (Cox 
et al. 2006). The present star-formation rate in 
NGC 3256 is estimated to be 33 AAq/yt from the 
far-infrared luminosity (Knierman et al. 2003). It 
is known that a tight relation between the most 
massive star cluster in a galaxy and the galaxy's 
star- format ion rate exists (e.g., Larsen 2004). In 
NGC 3256 we find about 10 clusters with masses 
in excess of lO^M©. If the star (cluster) forma- 
tion rate does increase substantially as the nu- 
clei merge, we may expect the NGC 3256 system 
to create clusters with masses significantly above 
10"^ Mq, such as those found in NGC 7252 and 
NGC 1316 (e.g., Schweizer & Seitzer 1998; Maras- 
ton et al. 2004; Bastian et al. 2006). 

Arp 220 is comparable to NGC 3256, as it 
also has an extremely high infrared luminosity, is 
thought to have formed through a merger, and 
has two distinct nuclei separated by only 300 pc 
in projection (Scoville et al. 1998, Wilson et al. 
2006). The nuclei in NGC 3256 are separated (in 
projection) by ^ 900 pc indicating that it is pos- 
sibly slightly younger (in terms of merger stage) 
than Arp 220. The star-formation rate in Arp 220 
(240 Mq/jt; Wilson et al. 2006) is approximately 
seven times higher than that in NGC 3256. It is 
then possible that NGC 3256 is currently poised to 
enter the regime of ultra-luminous infrared galax- 
ies (ULIRGs) as its star-formation rate increases 
due to the merging of the two nuclei. Arp 220 also 
closely follows the relation between global star- 
formation rate in the galaxy and the magnitude 
of the most massive cluster (Wilson et al. 2006), 
arguing further that NGC 3256 is going to form 
clusters in excess of lO^M©. 



5.3. Metallicities 

As Figure [9] shows the young clusters in the 
NGC 3256 system appear to have rather high 
metallicities, with the average being ~ 1.5Zq. 
This was also found for the clusters in the tidal 
tails described in Paper I. We do not see any 
major spread in metallicities among our clusters 
and, specifically, no differences coming from ages 
or emission versus absorption. The fact that the 
majority of these young clusters formed in a disk 
and their age spread is small fits in nicely with the 
notion of a normal starbust process. 

6. Comparison of NGC 3256 with Other 
Merging Galaxies 

Many of the clusters observed to be associ- 
ated with the molecular gas are quite massive 
(> lO^A^o), have survived for many internal 
crossing times (icr ~ 2-4 Myr), and are there- 
fore gravitationally bound. This justifies calling 
them young globular clusters. Such young glob- 
ulars have been found in many merging galax- 
ies, from beginning mergers (e.g., NGC 4038/39: 
Whitmore & Schweizer 1995) to completed merg- 
ers (e.g., NGC 1316: Goudfrooij et al. 2001a). The 
formation of these clusters is thought to trace the 
major star-formation events in these galaxies and 
they must form with approximately the same kine- 
matics as the gas out of which they form. It is 
therefore interesting to compare the NGC 3256 
cluster population to those of younger (in terms 
of dynamical stage) and older merging systems. 

The majority of the cluster population of the 
beginning merger NGC 4038/39 is still unam- 
biguously confined to the disks of the two pro- 
genitor galaxies (Whitmore et al. 2005; Bastian 
et al. 2006; Trancho et al. 2007). Therefore, 
NGC 3256 appears to predominately fall into this 
category (see § . 

Older systems, on the other hand, such as 
NGC 3921 (Schweizer, Seitzer, & Brodie 2004), 
NGC 7252 (Schweizer & Seitzer 1998), and 
NGC 1316 (Goudfrooij et al. 2001b) are all char- 
acterized by the kinematics of their clusters be- 
ing dominated by non-circular, halo-type orbits. 
When does the transition happen? It will be in- 
teresting to determine whether the majority of 
the star formation happens in the disks of the 
progenitors and their orbits are subsequently ran- 
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domized (i.e., turned into pressure supported sys- 
tems rather than rotational supported systems), 
or whether the star-formation events which pre- 
cede the destruction of the galactic disks pale in 
comparison with the star-formation rate during 
and after the destruction. 

Note that shortly after the merger a gaseous 
disk can reform around the nucleus of the merger 
remnant, which can harbor subsequent star for- 
mation (e.g., NGC 7252: Miller ct al. 1997; Wang, 
Schweizer, & Scoville 1992), although such star 
formation appears to occur at a much lower in- 
tensity than previous star-forming episodes during 
the merger. 

7. Summary and Conclusions 

We have studied the ages, metallicities, masses, 

extinctions, and velocities of 23 clusters in NGC 3256 
based on the Lick index system in conjunction 
with CO and HI maps. The main results are: 

• The clusters have rather high metallicities, 

with the average being ^ 1.5Zq (Fig. 9) and 
are massive, with masses in the range (2- 
40) xlO^ Mq. The ages of the clusters are 
between a few Myr and ~ 150 Myr. 

• There is strong evidence for a rotating 
molecular-gas disk in NGC 3256 (Sakamoto 
et al. 2006). The majority of the clusters in 
our sample follow the same rotation curve as 
the gas and hence were presumably formed 
in the molecular-gas disk. However, a west- 
ern subsamplc of five clusters has velocities 
that deviate significantly from the gas rota- 
tion curve. These clusters may either belong 
to the second spiral galaxy of the merger or 
may have formed in tidal-tail gas falling back 
into the system. 

• Although the merger began ~ 500 Myr ago 
(English et al. 2003), we found the clusters 
to be < 150 Myr old. Since there are still 
two distinct nuclei marking the presence of 
two galaxies, we conclude that the gas disk 
probably belongs to one of the galaxies and 
is not yet a disk of pooled gas produced in 
the merger itself. Presumably clusters older 
than the ones present in our sample do exist 
in NGC 3256. However, these older clus- 
ters would not have been selected as spec- 



troscopic candidates due their fainter mag- 
nitudes (i.e., only the brightest candidates 
were selected). 

• By comparing of the NGC 3256 cluster pop- 
ulation with other known galactic merg- 
ers, we suggest that this system is akin 
to Arp 220, although slightly dynamically 
younger. If this is the case, the wc may 
expect the star /cluster formation rate to in- 
crease significantly as the two galactic nuclei 
merge. This in turn may push NGC 3256 
into the category of ULIRGs (it is currently 
a LIRG). Due to the expected large increase 
the in the star/cluster formation rate, a few 
clusters above 10*" A40are predicted to form 
before this merger is through. 

• Some of the clusters which appeared as iso- 
lated in our ground-based images are clearly 
resolved into multiple sub-components in the 
HST-ACS images. The same effect has been 
observed in the Antennae galaxies, showing 
that clusters are often not formed in iso- 
lation, but instead tend to form in larger 
groups or cluster complexes. 

With these new results, i.e. cluster ages, 
metallicities, extinctions and kinematics, as 
well as recent CO and HI maps, N-body sim- 
ulations of this merger would be the best way 
to fully understand this wealth of data. The 
models would have important implications 
for (globular) cluster formation and destruc- 
tion as well as the star-formation history 
of the merger (through the age/metallicity 
of the clusters) with respect to other merg- 
ers hke the Antennae, NGC 7252 and Arp 
220. Finally, the details of the models may 
present important implications of the forma- 
tion of ellipticals galaxies through the major 
mergers of spiral galaxies. 
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Table 1 

Cluster properties. (The magnitudes have been only corrected for Galactic extinction) 
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^ 0=absorption, l=cmission, 2=cmission with WR features 
From Base position RA=10:27:00 DEC=-43:54:00 (J2000) 

These extinction have been calculated on spectra already corrected by Galactic extinction 
{Ay = 0.403) 

^ Ages calculated using the Starburst 99 models 
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Table 2 

Measured indices for the absorption line clusters. 
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Index definition by Scliweizer & Seitzer (1998). 
^ Lick index. 



Fig. 1.— ACS F555W image of NGC 3256 (cen- 
tral region) with the observed candidate cluster ID 
numbers. Green labels denote clusters with spec- 
tra dominated by emission lines, while red labels 
denote absorption line dominated cluster spectra. 
The line with an arrow points north, while the line 
without one points east. 
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Fig. 2. — Top: Histograms of ages derived by sim- 
ulating the effect of errors on the age-fitting rou- 
tine. The derived age and error (in logarithmic 
units) is given in each panel. Bottom: Same as 
top, except now for metallicity. See text for de- 
tails of the simulations. 
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Table 3 

Kinematics and Masses of the clusters. 
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Fig. 3. — Tests showing the effect of the S/N ra- 
tio on the measured hne indices. We used a high- 
S/N cluster spectrum, degraded its S/N ratio, and 
then measured line indices. Open symbols mark 
measurements made from the observed spectrum 
(corrected for extinction) directly, while solid sym- 
bols represent measurements made from a tem- 
plate spectrum derived for the cluster using the 
pPxF technique. Solid lines represent averages of 
the highest S /N experiments on the observed clus- 
ter spectrum (Clusterobs.oxt)- For further details, 
see text. 
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Fig. 5. — Examples of spectra for two clusters 
in our sample. The observed spectra have been 
corrected for the estimated interstellar extinction. 
The red lines represent the best fitting (see § 13.11 
for a discussion of the method) model template 
(age and metallicity). The green lines represent 
the residual (observed cluster — best fitting tem- 
plate — constant). The parameters for the best 
fitting template are given in each panel. 



Gonzalez-Delgado, gOOS: Geneva SSP Bruzual & Chariot (3003) SSP 




Fig. 4. — Determination of cluster ages and metal- 
licities. (left) H7 vs. [MgFe] from the Gonzalez- 
Delgado et al. (2005) SSP models for four differ- 
ent metallicities are shown. Data points with error 
bars mark observed clusters and their 1-cr errors. 
In addition, we show the position of the massive 
cluster W3 in NGC 7252 and the three tidal tail 
cluster in NGC 3256 from previous studies, (right) 
Same, but for Bruzual-Charlot (2003) models. 
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Fig. 6. — Example spectrum of an emission-line 
cluster, T2005, which also shows strong Wolf- 
Rayet features. 
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Fig. 9. — Metallicity distribution of the clusters 
in NGC3256. The figure includes the metallicities 
derived from both the absorption and emission line 
clusters. Note that all the clusters are fairly metal 
rich, with a mean around 1.5 Zq 



Fig. 7. — The position of the clusters in NGC 
3256 for different velocities (top left), ages(top 
right, where the ages are given in logarithmic units 
in years), extinctions (bottom left), and metallic- 
ities (bottom right). The contours are shown to 
highlight the main features of the galaxy. The 
upper and lower (magenta) circles mark the nu- 
cleus and the second brightest optical source in the 
galaxy, respectively. Triangles (red) mark clusters 
whose spectra are dominated by absorption lines 
while blue marks emission line clusters (sec also 
Table [1]). 
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Fig. 10. — HI position-velocity plot by English et 
al. (2003), with the observed clusters superposed. 
Open red circles mark the clusters that kinemati- 
cally follow the rotating CO disk, while filled blue 
circles mark those that do not follow the CO disk, 
but follow the HI velocities instead. 



Fig. 8. — A position- velocity diagram for the ob- 
served clusters in NGC 3256. The dotted line 
shows the rotation curve of the molecular g 
measured by Sakamoto et al. (2006). 
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A. Appendix A: Reduction of GMOS Data 

The data were reduced using the Gemini IRAF package. The individual tasks used at each of the reduction 
steps are listed below in parentheses. The raw GMOS images (3 CCDs read out using one amplifier each) are 
multi-extension FITS (MEF) files with one primary header unit, containing all the usual header information, 
and three pixel extensions, one for each of the detectors. 

A.l. Imaging Data 

Mean bias frames were created combining all the available bias frames taken during each of the GMOS-S 
observing runs (task gbias). 

Twilight flat fields were created from all the available twilight frames taken during each of the GMOS-S 
observing runs (gif lat). The overscan section of the images was trimmed off, and the images were bias- 
subtracted and flat-fielded using the task gireduce. The three image extensions in ciach exposure were then 
mosaiced into a single image extension in which the shifts and rotations between the three CCDs have been 
removed (gmosaic). 

The images were registered and co-added using the task imcoadd, and they have been calibrated using 
the science exposures taken under photometric conditions. 

A. 2. Spectroscopic Data 

We are using the MOS mode with approximately 25 spectra per mask. Each set of data include an 
arc and flat taken along with each science exposure and spectral dithers. Bryan Miller developed our own 
script-driven pipeline MOSPROC, based on Gemini IRAF scripts and customs IDL routines. The Gemini 
package scripts have been modified, mainly to improve the propagation of data quality planes. 

• MOSPROC 

— Bias substraction for spectra, fiats and arcs 

The bias subtraction for the spectroscopic data was done in the same way as for the imaging. 
Quartz-halogen flat flelds exposures were taken during the night, either after or before the science 
exposures. To run through MOSPROC, the science, arcs, and flats frames were first bias subtracted 
using the task gsreduce. 

- Bad Pixel Mask 

A bad pixel mask was constructed from the quartz-halogen fiats, taking into account the known 
bad pixels in each CCD, and masking the existing emission lines in the spectral QH lamps that 
came from the IR diffuser. 

- Flat-Field Correction 

The quartz-halogen flats were overscan trimmed and bias-subtracted like the imaging data and 
the mask definition file (MDF) was added as a table extension to the MEF. The MEF contains 
the locations of the slits in the focal plane and is used for bookkeeping during the remaining 
reduction. The flats were normalized by fitting a high order polynomial to each line to remove 
the shape of the quartz lamp but leaving the fringe pattern at the red end of the spectra gsf lat. 
The task gsreduce was then used to divide the science data by the flat fields. 
The resulting science exposures still have three extensions, but now trimmed, bias-subtracted, 
and flat-fielded. 

— Wavelength Calibration and Distortion Correction 

The wavelength calibration was determined from CuAr lamp spectra taken either before or af- 
ter the science exposures. The dispersion function was fit with a sixth-order polynomial that 
gave a typical rms error of 0.3 Angstroms (gswave length). Within each slitlet the position of 
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each arc line with spatial position is used to rectify and wavelength calibrate each 2D spectrum 
(gstransf orm). 

Cosmic Rays Removal 



Cosm ic rays were identified and removed using a Laplacian edge detection algorithm (jvan Dokkum 



20011 ). The locations of the cosmic rays are recorded in the 2D data quality image for each slitlet. 
The spectra were then traced, background-subtracted, and extracted using gsextract, a wrapper 
for twodspec . apall, to allow the handling of MEF files. The spectra were traced using a 5 order 
polynomial and averaging every 50 pixels in the dispersion direction. Background subtraction 
was done by fitting a second-order polynomial with 3<t rejection to a region perpendicular to the 
trace. The variance and background spectra were saved to help with the error estimation. 

— Quantum Efficiency (QE) Correction 

The quantum efficiency (QE) as a function of wavelength for the three CMOS CCDs can differ by 
up to about five percent at the given wavelength. If this difference is not corrected then spectra 
can have noticable jumps at the gaps between the CCDs. We have measured the relative QE of 
CCDl and CCDS with wavelengh compared with CCD2 using QH flats taken at 25nm intervals 
from 350nm to 700nm. The corrections are applied to the data using an IRAF script called 
qecorr. 

Correction for Slit Losses 

The CMOS instruments do not have atmospheric dispersion correctors and in general the MOS spec- 
tra were not taken with the slits parallel to the parallactic angle. Therefore, there are wavelength- 
dependent slit losses due to the difference between the parallactic angle of each exposure and the 
position angle of the slits. In order to combine the spectra from different exposures properly these 
difference must be corrected. We calculate the slit losses based on measu red image quality, slit width 
(0.75 arcsec), and the PA-parallactic angle difference using the method of iFilippenk 



Relative Flux Calibration and Reddening 

— Relative Flux Calibration 

The relative instrument spectral response function was determined using observations of a flux 
standard star. Spectra were taken at central wavelengths of 400, 500, and 600nm so that the 
combined sensitivity function covered the full wavelength range of the MOS spectra. The standard 
spectra were reduced in exactly the same way as the MOS spectra. The sensitivity function was 
computed using the IRAF task gssensfunc and this was applied to the MOS spectra using the 
task gscalibrate. 

As the last step, we combined the individual spectra of the same sources using our IRAF task 
called gscombine. This uses the final data quality planes to mask the chip gaps and other bad 
pixels and then scales and averages the spectra using scombine. 

— Reddening 

The interstellar extinction along the line of sight towards NGC 3256 is Ay — 0.403. The final cali- 
brated spectra were corrected using the empirical selective extinction function from lCardelli et al.l 
included in the task noao . onedspec . deredden. 



G.T. would like to thank to Matt Mountain and Phil Puxley for the tremendous support throughout this 
project. 
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